Quantum computing with spin cluster qubits 
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We study the low energy states of finite spin chains with isotropic (Heisenberg) and anisotropic 
(XY and Ising-like) antiferromagnetic exchange interaction with uniform and non-uniform coupling 
constants. We show that for an odd number of sites a spin cluster qubit can be defined in terms 
of the ground state doublet. This qubit is remarkably insensitive to the placement and coupling 
anisotropy of spins within the cluster. One- and two-qubit quantum gates can be generated by 
magnetic fields and inter-cluster exchange, and leakage during quantum gate operation is small. 
Spin cluster qubits inherit the long decoherence times and short gate operation times of single spins. 
Control of single spins is hence not necessary for the realization of universal quantum gates. 
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Quantum computers outperform classical computers 
on certain tasks [1-4]. The main challenge on the way to 
a universal quantum computer is to achieve control over 
single quantum mechanical two state systems (qubits) 
while preserving long decoherence times. Electron [5,6] 
and nuclear [7,8] spins have been identified as promis- 
ing candidates for qubits because they are natural two 
state systems and decoherence times for the spin degree 
of freedom are unusually large [9,10]. 

For both electron [5] and nuclear spin [8] qubits, one- 
qubit gates can be realized by local magnetic fields or 
by electrically tuning a single spin into resonance with 
an oscillating field. Two-qubit gates rely on electrical 
control of the exchange interaction between neighboring 
electron spins. However, even for electrons in quantum 
dots with a typical diameter of 50 nm, the required local 
control over electrical and magnetic fields is challenging. 
One possibility to circumvent the problem of either local 
magnetic fields [11] or local exchange interaction [12] is 
to encode the qubit in several spins. Such encoding has 
also been studied in the context of coherence-preserving 
qubits [13]. However, all these schemes still require con- 
trol at the single-spin level. 

More generally, the requirements on both local mag- 
netic and electrical fields can be relaxed by increasing 
the size of the qubit. In the present work, we show that, 
for a wide class of antiferromagnetic spin s — 1/2 chains 
with an odd number of sites, n c , 
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the ground state doublet of the array [Fig. 1(a)] can de- 
fine a new "spin cluster qubit" for which quantum gate 
operation times and decoherence rates increase only mod- 
erately with array size. These features are surprisingly 



stable with respect to anisotropy (Jj_ ^ Jz) and spatial 
variation (described by ft) of the intracluster exchange, 
the spatial shape of the fields controlling quantum gate 
operation, and the cluster dimension. Spin cluster qubits 
can be realized in a wide variety of systems, e.g. ar- 
rays of quantum dots [5,14], clusters of P atoms in a 
Si matrix [8], and electron spins in molecular magnets. 
In contrast to the encoded qubits suggested in earlier 
work [11-13], quantum computation with spin cluster 
qubits is possible without control over local spin inter- 
actions. 

Isotropic spin chains as qubits. - For illustration, we 
first discuss a spin chain with isotropic uniform exchange, 
J± = J z > and fi = 1 in Eq. (1). Energy eigenstates 
can be labeled according to their quantum numbers of to- 
tal spin S = 53i=i ^* an d f ne ^-component of total spin, 
S z . Because of the antiferromagnetic exchange, states in 
which the total spin of the chain is minimized are energet- 
ically most favorable [15]. For odd n c , the ground state is 
a S = 1/2 doublet separated from the next excited state 
by a gap A ~ J7r 2 /2n c determined by the lower bound 
of the des Cloiseaux-Pearson spectrum. We define the 
spin cluster qubit in terms of the S = 1/2 ground state 
doublet by S z \0) = (h/2)\0) and S z \l) = -(h/2)\l). The 
states {|0),|1)} do not, in general, have a simple rep- 
resentation in the single spin product basis, but rather 
are superpositions of n c !/[(n c — l)/2]![(n c + l)/2]! states 
[Fig. 1(b)]. For example, for the simplest nontrivial spin 
cluster qubit with n c = 3, 
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and |1) is obtained by flipping all spins. 

In spite of their complicated representation in the 
single-spin product basis, |0) and |1) are in many respects 
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very similar to the states | f) and | J.) of a single spin and, 
hence, can be used as qubit states for universal quantum 
computing [5]. Because {|0), |1)} belong to one S = 1/2 
doublet such that 5"|0) = h\l), and = fi|0) where 

= S x ± iS y , a magnetic field B constant over the 
cluster acts on the spin cluster qubit in the same way 
as on a single-spin qubit. Hence, both the one-qubit 
phase shift and the one-qubit rotation gate can be gen- 
erated by magnetic fields B z (t) and B x (t), respectively, 
possibly in combination with (/-factor engineering [5,16]. 
For a given B ZiX (t), operation times of one-qubit gates 
are equal to the ones for the single-spin qubit. We note 
that, due to quantum mechanical selection rules, we have 
(i|S|0) = (i|S|l) = for \i) ^ |0),|1), i.e., a uniform 
magnetic field does not cause leakage to states outside 
the ground state doublet. 

For the CNOT gate, one requires a tunable exchange 
interaction between one or several spins of neighbor- 
ing spin cluster qubits I and II. For simplicity, we first 
restrict our attention to an isotropic exchange coupling 
H* = J»(t)s^ c • between the outermost spins of clus- 
ters I and II, respectively [Fig. 1(c)]. This exchange in- 
teraction will in general not only couple states within the 
two-qubit basis {|00), |01), |10), 1 11)}, but will also lead 
to transitions to excited states (leakage) . If J* (t) changes 
adiabatically, i.e., on time scales long compared to ft/ A 
and |J*(i)| -C A for all times t, leakage remains small 
(see below). The action of can then be described by 
an effective Hamiltonian in the two-qubit product basis 

H* = J* z {t)S l z S l z l + ^M(S I+ S U - + S'-S 1 ^), (3) 

where the roman numbers label the spin clusters, 
J„(t) - 4J,(t)| I (0|4 c , z |0)i||n(0|J I 1 i z |0)„|, and J. ± (i) = 
4J.(t)|i<l|^ e ^0)i||n<0|«£ x |l)ii|. We have shown that 
the coupling is isotropic also in the two-qubit prod- 
uct basis and acts on the states |0) and |1) of neighboring 
spin chains in the same way as an isotropic exchange be- 
tween two single spins. |i(1|sJj c;X |0)i| and |ii(0|s" x |1)ii| 
determine the gate operation time tcn of the CNOT gate, 
\i) \j) — > \i) \i+j mod 2) where i, j = 0, 1. For n c = 9, . . . , 
15, the matrix elements are of order 0.1, i.e., a factor of 
5 smaller than for a single spin 1/2. 

Although we have so far discussed one-qubit gate 
operations induced from spatially uniform magnetic 
fields, such uniformity may be difficult to achieve 
experimentally. One-qubit gates can be performed 
with spatially varying fields B ix and B iz (or g fac- 
tors) for which |(1| Yh=i 9i^BB itX s itX \0)\ ^ and 
(0| Y17=i Qif^BBi z §i z \0) ^ 0, respectively. Similarly, the 
analysis leading to Eq. (3) remains valid for a wide class 
of coupling Hamiltonians H* for which (10 1 H* | 01) ^ 0. 
For illustration we discuss two examples. First, couplings 
between several spins of cluster I to several spins of clus- 
ter II, such as = J* Y^i=i si • sj 1 , are permitted and 



even lead to a decrease of tcn because the coupling of 
several spins in the microscopic Hamiltonian leads to an 
increased effective coupling between the clusters. Second, 
a modification of the intracluster exchange couplings by 
H* due to additional terms such as J^s^ 1 • does not 
invalidate the proposed gate operation scheme. This il- 
lustrates the most significant advantage of the spin clus- 
ter qubits over single-spin qubits - that it is sufficient 
to control magnetic fields and exchange interactions on 
a scale of the spin cluster diameter. For the linear spin 
cluster qubit, this length scale is n c times larger than the 
original qubit. 

A set of universal quantum gates is necessary but 
not sufficient for the realization of a quantum computer. 
Rather, additional requirements must be met, including 
initialization, decoherence times large compared to gate 
operation times, and readout [17]. Initialization can be 
achieved by cooling in a magnetic field B z to a temper- 
ature [5] T < gnBB z /kB < A/fc^. Because the state 
of the spin cluster qubit, |0) or |1), determines the sign 
of the local magnetization at each site within the spin 
chain [Fig. 1(b)], readout of the spin cluster qubit can 
be accomplished by readout of the spins within the clus- 
ter [5,18]. 

An important consideration is the effect of decoherence 
on spin cluster qubits. The scaling of the decoherence 
rate with system size depends on the microscopic deco- 
herence mechanism. For electron spins in quantum dots, 
fluctuating fields and nuclear spins have been identified 
as dominant sources [5,6,19]. We model [5] the action of 
fluctuating magnetic fields by = b(t)S z where b(t) 
is Gaussian white noise, (b(t)b(0)) = 2Tr-f B S(t). Because 
the magnetic moment ±g/zs/2 of the spin cluster qubit 
is the same as for a single spin, the decoherence rate [20] 
7T7 S is independent of n c . In contrast, the decoherence 
rate due to fluctuating fields acting independently on 
each site increases linearly with n c . 

Spin dynamics during gate operation. - One- and two- 
qubit gates induce spin dynamics in the clusters, and 
leakage out of the ground state doublet is required to 
remain small. In order to quantify leakage, by numerical 
integration of the Schrodinger equation we trace the time 
evolution of a small spin cluster qubit (n c = 5) during the 
one-qubit rotation gate. The qubit is rotated coherently 
from |0) into |1), which corresponds to a simultaneous 
rotation of all spins [Figs. 2(a) and (b)]. The one-qubit 
rotation can also be generated by an inhomogencous field 
B x acting, e.g, only on the central spin of the cluster as 
long as g/j,BB x <C A [Figs. 2(a) and (b)]. Leakage due to 
instantaneous switching is less than 0.3% for g/iBB x = 
0.1J, but increases with g\xsB x [Fig. 2(b)]. 

For the special cases J* z = J*± and J %z = in 
Eq. (3), an explicit pulse sequence for the CNOT gate 
has been derived previously in Refs. [5,21]. We de- 
fine the unitary time evolution operator C/*(7r/2) = 
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Texp(-ifdtH*/h), with - / dtJ*±(t)/h = tt/2. 
Then, more generally, 

U C NOT ~ e- l7rS "/2 e ^n 1 .SV3 e42 ™ 2 .S"/ 3;74i(7r/2) 

x e™ s l U* (ir/2)e-™ s */ 2 e^ sl ^ 2 e™ s y/ 2 (4) 

is the CNOT gate for an arbitrary effective XXZ 
coupling Hamiltonian [Eq. (3)] if J*j_ 7^ 0, where ni = 
(l,-l,l)/\/3 and n 2 = (1,1,-1)/V3 . We confirmed 
that for the complete pulse sequence the dynamics of 
two spin clusters is as predicted on the basis of the two- 
level description [Fig. 3] . Leakage induced by is small 
for </* <C A because all spins in the clusters corotate 
although couples only the outermost spins. 

Spatially varying and anisotropic exchange. - We 
show next that spin cluster qubits are extremely robust 
with respect to spatial variation (accounted for by fi in 
Eq. (1)) and anisotropies (Jj_ 7^ J z ) of the intraclus- 
tcr exchange. For spatially varying isotropic exchange 
(J_i_ = J z = J), the system still exhibits a S = 1/2 
ground state doublet [15] and the above analysis remains 
valid. In systems such as quantum dot arrays where it is 
possible to engineer the intracluster exchange J fi during 
sample growth, the qubit basis states {|0),|1)} can be 
tailored to some extent. 

We next consider the XY chain, J z = 0. By the 
Jordan- Wigncr transformation [22], the XY spin chain 
is mapped onto a system of noninteracting spinless 
fermions with spatially varying hopping amplitudes, H = 
-(^-l/ 2 ) Ya=i 1 fi(i>i+ii>i where tpi annihi- 

lates a Jordan- Wigner fermion at site i. We find that the 
one-particle Hamiltonian has (n c — 1)/2 states with nega- 
tive and positive energy, respectively, which are pairwise 
related to each other by staggering of the wave function. 
There is one zero-energy eigenstate 

e cx 1,0,- — ,0, — — ,0,...,±— . (5) 

\ ,/2 J2J4 j2/4 • ■ • Jn c -1 J 

The ground state doublet of the XY chain corresponds 
to the lowest (n c — l)/2 and (n c + l)/2 Jordan- Wigncr 
fermion levels filled. For fi = 1, A ~ wJ±/n c . Similarly 
to the spin chain with isotropic exchange, one-qubit gates 
can be realized by magnetic fields B z (t) and B x (t) unless 
(1|&|0) = 0. For n c < 9 and f t = 1, |(l|£ x |0)| > 0.4. 
From Eq. (5), one can also calculate all matrix elements 
entering Eq. (3). In particular, for fi = 1, (0|s„ cjZ |0) = 
l/(n c + l), and |(l|sn c ,x|0)| = l/y/2(n c + 1). Because 
of the anisotropy of the intrachain exchange, -ff* (which 
is isotropic in the single spin operators) translates into 
an anisotropic effective Hamiltonian Eq. (3). Never- 
theless, the CNOT gate can still be realized accord- 
ing to Eq. (4). For the anisotropic chain, a magnetic 
field applied along an axis n translates into a rota- 
tion around the axis cx (2| (llS^O)!?^, 2| (l|»S , 2: |0)|n y , n z ) 



in the Hilbcrt space spanned by {|0),|1)}. A one 
qubit rotation around an arbitrary axis hence requires 
appropriate rescaling of B. For example, the rota- 
tion corresponding to exp(i27rni • S ! /3) [Eq. (4)] for the 
isotropic chain can be achieved by applying a magnetic 
field B = B Q (1 + 2/(21(115*10) |) 2 ) 1/2 /\/t along the axis 
cx (1,-1,2|(1|S X |0)|) for a time 2irh/3gLi B B . For given 
J* and B, the CNOT gate operation time increases at 
most linearly with n c . 

For J z » J± (Ising-like systems), where 
|0) = I T)i| D2 • • • I t)„ e +0{J±/J z ), the ground state 
doublet is separated from the next excited state by 
an n c -indepcndcnt A ~ J z min(fi). In perturbation 
theory in J±/J z , for fi = 1, the matrix elements 
|(l|^x|0)|,|(l|a nc ,x|0)| ~ (2J ± /J,)(^-i)/2 decrease ex- 
ponentially with n c . Even for medium sized chains n c > 9 
and J±/J z < 0.2, an isotropic translates into an ef- 
fective Ising Hamiltonian, J*± ~ in Eq. (3). Hence, 
only quantum computing schemes which rely on Ising in- 
teractions [23] are feasible. 

Discussion. - The main idea of the present work ap- 
plies not only to spin chains but remains valid for a wide 
class of antiferromagnetic systems with uncompensated 
sublattices, also in higher dimensions d > 1 and for larger 
spins s > 1/2. Wc illustrate the advantages of spin clus- 
ter qubits for electron spins in quantum dots with a typ- 
ical diameter of d = 50 nm, where the exchange cou- 
pling can be as large as 10 K [6]. One-qubit operations 
are realized, e.g., by (/-factor engineering in presence of 
a static field B ~ 1 T. We now compare the perfor- 
mance of a spin cluster qubit formed by n c = 5 quan- 
tum dots coupled by an intracluster exchange J = 10 K 
with a single spin qubit. To obtain an estimate, we con- 
sider gate operation by switching the magnetic field B 
to g^BB = 0.7 K, and J* = 2.3 K [6], small compared 
to the energy gap A = 7.2 K of the spin cluster. For 
single spins, the gate operation times for the NOT and 
CNOT gate are 36 ps and 117 ps, respectively. Assum- 
ing that the magnetic field decreases smoothly from its 
maximum value at the central spin of the spin cluster 
qubit to 0.2_B acting on spins 2 and 4, we find that the 
operation time for one-qubit gates increases by a fac- 
tor l/2\{l\s 3tX + 0.2s 2 , x + 0.2s 4 , x |0)| = 2.2 compared to 
the single spin. Similarly, for the operation time of the 
CNOT gate we find 280 ps for the spin cluster qubit. The 
main advantage of the spin cluster qubit is that it is suf- 
ficient to control magnetic fields or g factors on a length 
scale of n c d = 250 nm and exchange couplings on a scale 
of 2n c d = 500 nm. This would allow one to control the 
exchange between neighboring clusters optically [14] at 
the expense of an increase in gate operation times by a 
factor of 2. 

Other possible applications for spin cluster qubits in- 
clude, e.g., P atoms in a Si matrix [8] and molecular 
magnetic systems [24] . For electron spins of P atoms in 
a Si matrix, the requirement of positioning P with lat- 
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tice spacing precision [8,25] can be circumvented by the 
use of spin clusters instead of single spins. More gener- 
ally, the present work shows that for universal quantum 
gates control is not required at the level of single electron 
spins. Because a qubit can always be mapped onto a spin 
1/2, the general principle of arranging several qubits into 
a cluster qubit applies to any quantum computing pro- 
posal. 
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FIG. 1. (a) The states of the spin cluster [Eq. (1)] define 
the spin cluster qubit. (b) |0) and |1) have a complicated 
representation in the single-spin product basis, as evidenced 
by the local spin density, (c) Quantum gates are generated by 
magnetic fields or gr-factor engineering (one-qubit gates) and 
a switchable inter-qubit exchange J*(t) (two-qubit gates). 
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FIG. 2. (a) Local spin density within a spin cluster qubit 
(n c = 5) as function of <j) oc gfiBB x t/h obtained by integration 
of the full Schrodinger equation for homogeneous (solid line) 
and inhomogeneous B x (dashed and dashed-dotted lines), (b) 
For B x A/g/j,B or homogeneous B x , the state rotates co- 
herently from |0) to |1). For a magnetic field acting only on 
the central spin of the cluster, the leakage increases to 7% for 
gfi B B x = 0.5J. 
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FIG. 3. CNOT gate for two small spin cluster qubits 
(n c — 3) obtained by numerical integration of the Schrodinger 
equation [see Fig. 1(c)]. The plotted probabilities and the 
phases (not displayed) show that (a) |00) — > |00) and (b) 
10) — * 1 11). We have chosen a pulse sequence [Eq. (4)] with 
instantaneous switching (at times ti), B = O.lJ/g^B , and 
J» = 0.1J. Leakage due to instantaneous switching (0.7% for 
our parameters) can be reduced by decreasing J* and B. 
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